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Abstract. It is shown that like Bianchi I, V and IX models, a Kantowski-Sachs 
cosmological model also allows a unitary evolution on quantization. It has also been 
shown that this unitarity is not at the expense of the anisotropy. Non-unitarity, 
if there is any, cannot escape notice in this as the evolution is studied against a 
properly oriented time parameter fixed by the evolution of the fluid. Furthermore, we 
have constructed a wave-packet by superposing different energy eigenstates, thereby 
establishing unitarity in a non-trivial way, which is a stronger result than an energy 
eigenstate trivially giving time independent probability density. For a ^ 1, we have 
proved that the Hamiltonian upon suitable operator ordering admits a self-adjoint 
extension under reasonable assumption, using a standard theorem. This consolidates 
the recently shown fact that the problem of non-unitarity cannot be a generic problem 
of anisotropy. 
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1. Introduction 

The Universe as a whole, being a viable physical system, is expected to have an un¬ 
derlying quantum theory. The interaction that governs the Universe at large scale is 
gravity, which as yet does not have any consistent and universally accepted quantum 
theory and thus one cannot really look at the quantum universe as an application of 
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quantum gravity. The alternative, an attempt to bypass the conceptual difficulties of 
quantum gravity, is provided by Wheeler-deWitt equation[l, 2, 3], which is actually the 
starting point for quantizing a cosmological model, and forms the subject of quantum 
cosmology. However, quantum cosmology by itself is not actually free from conceptual 
problems as well as some practical difficulties. One raging problem is that of the choice 
of a properly oriented time parameter, as in a relativistic theory time is a coordinate 
and not a scalar parameter same for all observers. This problem, and the possible res¬ 
olutions in various ways constitute almost a branch of quantum cosmology [4, 5, 6, 7]. 
Recently, a choice of an oriented time parameter depending on the evolution of a per¬ 
fect fluid constituting the matter distribution in the universe has been in use with a fair 
amount of success. The strategy actually had been given long back by Lapchinskii and 
Rubakov[8] based on Schutz’s formalism of writing the velocity vector of the fluid in 
terms of some thermodynamic potentials[9, 10]. Another serious conceptual problem, 
that of the interpretation of the wave function and that of the boundary conditions are 
reviewed by Wiltshire [11], Halliwell[12] and very recently by Pinto-Neto and Fabris[13]. 

One nagging problem of quantum cosmology had been the fact that anisotropic 
models are found to be nonunitary and hence do not respect the conservation of proba¬ 
bility. As the collection of all the 3-space Riemannian metric, the Superspace, is too big 
to work in, the usual practice is to restrict the description in a truncated version of that, 
namely a minisuperspace. The natural choice of the minisuperspace is indeed a spa¬ 
tially homogeneous and isotropic FRW 3-space which is compatible with the presently 
observed universe. But anisotropy has its relevance. The formation of the present 
structure of the Universe requires a small but hnite anisotropy in the CMB tempera¬ 
ture ~ 10“^) and the recent observations also reveal the same order of magnitude 
of anisotropy. The other requirement is aesthetic; there is no reason, a priori, for the 
universe to be isotropic. So a correct version of a quantum theory should not lead to 
inconsistent physics, such as a nonconservation of probability, even for the anisotropic 
models. It is interesting to note that without the existence of a time parameter with the 
correct orientation, this nonconservation of probability in quantum anisotropic models 
is generally obscure and may escape being detected[14, 15]. This alleged pathology of 
non-unitarity often ascribed to the hyperbolicity of the Hamiltonian of the anisotropic 
cosmologies [16]. 

Very recently it has been shown that the nagging problem of non-unitarity can 
actually be cured. For a Bianchi type I model, either by a clever ordering of opera¬ 
tors or by a suitable transformation of coordinates, one can hnd self-adjoint extension 
of the Hamiltonian[17]. In fact this is not exclusive for a Bianchi I model where the 
spatial curvature is zero, the similar self-adjoint extension is possible for anisotropic 
spaces with negative and positive spatial curvature given by Bianchi V and IX models 
as well[18]. It deserves mention that both these investigations involve a properly ori¬ 
ented time parameter through the evolution of the fluid as mentioned earlier. Of course 


Unitary evolution for a quantum Kantowski-Sachs cosmology 


3 


the construction of the wave packet and the evaluation of the norm depends on the 
degree of difficulty of the integration involved, so one cannot work out what may be 
called the “general” situation, one actually has to depend on some specihc equations of 
state for the fluid present. But one counter example is good enough to disprove the fact 
that alleged non-unitarity is a generic problem of anisotropic models, and one now has 
several nontrivial examples. It deserves mention that even istropic models actually can 
show a nonunitary evolution unless the operator oredring is carefully chosen. Issues of 
unitarity of a quantized isotropic model in the presence of a scalar held has been very 
recently discussed by Almeida etal[19]. 

The motivation of the present work is to check if unitary evolution is a possibility 
in a Kantowski-Sachs (KS) cosmology, which is qualitatively different from Bianchi I, 
V and IX models. The KS metric represents a homogeneous but anisotropic spacetime 
and has its own relevance in gravity in various ways, such as the fact that a KS metric 
is isometric to the metric for the interior of a spherically symmetric black hole. For a 
concise list of the relevance of a KS metric, we refer to the recent work by Parisi, Radi- 
cella and Vilasi[20]. A quantization scheme for a KS metric, although not much talked 
about, is not completely new. It has been discussed in the context of a non-commutative 
geometry by Garcia-Compean, Obregon and Ramirez[21], and very recently in the con¬ 
text of a loop quantum cosmology by Joe and Singh[22]. Quantizing a KS cosmology in 
a more standard version of gravity was discussed long back by Conradi[23]. The latter 
includes a fluid, namely a pressure-less fluid in the quantization scheme, but does not 
indicate anything regarding unitarity. 

We show that the Hamiltonian is indeed self-adjoint in the case of a Kantowski- 
Sachs model with a stiff fluid (i.e. fluid with an equation of state P= p, where P is 
pressure and p is density of fluid), upon choosing a suitable weight factor while dehn- 
ing the norm. An explicit example of a unitary solution has also been obtained. The 
complete calculations could be possible, however, only for a particular equation of state, 
that of a stiff fluid. For other equations of state, we can not hnd any explicit solution 
owing to the computational difficulty. Nonetheless, we could show that even for such 
cases {P 7^ p)the Hamiltonian admits a self-adjoint extension under reasonable assump¬ 
tion with a clever ordering of operators. 

The general formalism and the Wheeler deWitt equation for the KS model is given 
in section 2. Section 3 takes up the case of the particular example with a stiff fluid. In 
section 4, we arrive at the relevant Wheeler deWitt equation for the general barotropic 
fluid and show how a suitable operator ordering can make Hamiltonian self-adjoint 
extendible. Section 5 includes some discussions. 
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We start with the standard Einstein-Hilbert action 

[ d'^xy/^R + 2 [ ^fhhabK^'^ + [ d‘^x-\/^P, (1) 

J M JdM Jm 

where is the extrinsic curvatnre, the induced metric over the boundary dM of 
the 4 dimensional space-time manifold M and P is the fluid pressure, given by P = ap, 
where p is the density of the fluid and a is a constant. The units are so chosen that 
levrG = 1. 

The Kantowski-Sachs metric is given by 

ds^ = —n^df^ + X‘^{t)dr‘^ + \d6‘^ -I- sin^ {9) , (2) 

where n{t) is a lapse function while X{t) and Y{t) are scale factors. In order to 
facilitate the process of quantization, we go over to the Misner representation via the 
transformation [24] 

X{t) = e^^+, ( 3 ) 

r(t) = (4) 

. This transformation allows us to write the metric (2) as 

ds^ = -n^de + + e -2 V 3(/3++/^-) [^^2 ^ 


The fluid sector can be expressed in terms of some thermodynamic variables following 
Schutz formahsm[9, 10], which has been subsequently developed by Lapchinskii and 
Rubakov[8] in context of quantum cosmology. The method has been used quite 
extensively later by many [17, 18, 25, 26, 27]. In this formalism, the Lagrangian for 
the system can be extracted out of the action (1) as 


n 


'/?+ - 

+2ne^^+ + 




n “ - 


a 


(l + «)' 






The metric as well as all other quantities are spatially homogeneous, the integration of 
space part yields a constant in (1) and is thus inconsequential as it can be absorbed 
in the right hand side zero of the variational principle. The quantities d, e, s are 
thermodynamic potentials which determine the velocity vector. The details are given 
in [8]. The last term within the square bracket is the contribution from the fluid 
sector whereas the rest denotes contribution from the gravity sector. The corresponding 
Hamiltonian for the gravity sector for the metric (5) can be written as 


H - ILpV3{P++2P-) 

" “ 24 


2 I 2 

-Pp_ + P/34 


-48e 


- 2 ^ 3 / 3 - 


where the canonical momenta are defined in the usual way. 
momenta for the fluid sector as and Ps = ^ 

Hamiltonian is written as 


( 7 ) 

We define the canonical 
and the corresponding 


Hf = 


( 8 ) 


cole 
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With the canonical transformation, 

T = -psexp{-S)pr-\ (9) 

PT = pf^^expiS), ( 10 ) 

e' = e + (a + 1) —, (11) 

Pe 

Pe=Pe, (12) 

the Hamiltonian for the fluid sector becomes 

Hf = (13) 


Hence combining (7) and (13), the super Hamiltonian is given by 


gV3(l-a)(/3++2/3_) r ^ ^2^ _ + 24 ^^ 114 ) 


H = Hg + Hf 

^g«x/3(/J++2/3_) 

“ 24 

Now, one can choose a gauge n = 24noe““'^^^+'''^^“^ and vary (14) with respect to 
uq to obtain the Hamiltonian constraint as 

gv/3(l-a)(/3++2/3_) | _p2^ _ 486"^^^- | + 24pr = 0. (15) 


3. Stiff fluid (p = P) 


As an example, we choose a stiff fluid given by a = 1. This choice avoids the ordering 
ambiguity since = 1 for a = 1. On quantizing the model, we have the 

Wheeler-deWitt equation. 


<92 




48e"^^^- I ?/> = 24^ 


dfj 

’In'’ 


(16) 


in units of h = 1. The canonical momenta Pj 3 ^ and pt are replaced by and 

respectively following the procedure elucidated in [16, 27, 17]. 


3.1. Probability density 

With the separability ansatz fj = e*'^^+^+^(/5_)e“*®^, (16) becomes 

{ ^ + 3^+ - 48e-2^^-1 0 = 24770, (17) 

which has the solution 

0 = 77,, (4e-^^-), (18) 

where u = \/kf — 877. It deserves mention that this solution is similar to that obtained 
by Garcia-Compean et a/[21]. The expression for u in the present work is different due 
to the contribution from the fluid sector. 
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Now, we note that sector admits a plane wave like solntion. Hence, we might be 
tempted to form a Gaussian wavepacket superposing different plane wave solutions by 
doing an integration over A:+ variable. But this is actually not that simple. If we want 
to construct a wavepacket with dehnite energy E, then varying would imply varying 
u = \Jk\ — 8E, hence we are forced to superpose j ’s as well, resulting to 

a complicated wavepacket, which is not Gaussian at all. This implies that if we intend 
to form a Gaussian wavepacket corresponding to /5+ sector, it can not be an energy 
eigenfunction, which motivates us to dehne a new parameter 

a = k"]^ — SE = iS' (19) 

and trade off (i?, k+) against (a, kjf) i.e the wavefunction can be thought of dependent 
on and a independently, and E becomes a parameter dependent on and a through 
equation (19). 

In terms of a, the wavefunction for (3- sector reads as 

( 4 %), ( 20 ) 

where we have dehned, for brevity, y = e~^^-. 


Now we can superpose the plane wave like solutions for /5+ sector , keeping a hxed. 
This, in turn, means the resulting wavepacket will no more be an energy eigenstate, it 
is rather a superposition of energy eigenstates. The wavepacket is given by following 
expression. 


= <^(x)e*«'^0(/5+), 

(21) 

where we dehne 0(0+) as 


0(0^) = J 

(22) 


If we want to insist on superposition of negative energy eigenstates only, the inte¬ 
gral over /c+ should be in the interval y^).The requirement of a negative E for 

the gravity sector stems from the demand that the energy for the fluid part is positive 
whereas the total energy is zero in view of the Hamiltonian constraint. Similar kind of 
consideration has been made in earlier work [17, 18]. 


Upon integration over the said interval, we have 


rf = . =\/7rexp 


\/—8i + T 
where g{a) is given by 


-/c^T 8\/3 A:+o/5+ + 

-8i + T 


g{a) = Erfi 


^ (8/c+o + V^{-a - iT) + 4^\/3/5+) 


\/—8z -|- T 

With the dehnition, 11-011 = we have 


[g{8)-g{-8)] (23) 


( 24 ) 
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m 


1 9 

—= / 

VsJ-^ 


1 fn 
~ 2 


Erf 


\fif\fa + /c+o) 


Erf 


— + /c_|_o) 


(25) 


We have specifically shown, for the (3+ sector, that the norm is hnite and time 
independent. This happens non-trivially since we are not dealing with energy eigenstate. 
We recall that for energy eigenstate, the probability density function is trivially time 
independent. Now we turn our attention to /3_ sector. We recast the equation (17) 
using the variables x and a\ 


^ TT + -a)(j) = 0 

dx dx 

This equation can be rewritten in the standard self-adjoint form: 


A. 

dx 


X 


d(j) 

dx 


+ 16x 


= 0 


with definition of inner product is given by 

POO 

{<t'i\4>2) = dxx <l>]{x)Mx) 


(26) 


( 27 ) 


(28) 


Hence, the Hamiltonian for /?_ sector is self-adjoint as well ensuring a unitary time 
evolution. The norm can explicitly be calculated and shown to be hnite, since 

/»rv^ 

TIU 


dx X (p*{x)(p{x) = 


32 sinh (ttu) 


(29) 


where (j){x) = (4%) and = a. 


The relevant plots of probability density function is given below at different time 
instant. We note that, f is time independent, hence the probability density due to [3- 
sector (we will call it Pr_ = X'f{x)f'*{x)) is time independent, while probability density 
due to /9+ sector, which we call = -^(/3+)?/)*(/3+), is time dependent (even though 
the norm is hnite and time independent). 



Probability Density at T=6, o-=2, *+(,= 1 Probability Density at T=12, <9=2, *+0=1 

Pr+ "+ 



Figure 1: Probability Density {Prjf) 

We note that the wavefunction dies out for large y, which translates to the fact that 
the probability of having (3- = —oo is zero. Similarly the probability density dies out as 
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Probability Density, (T=2, fixed 13+ 
Pr_ 



Figure 2: Probability Density (Pr_) 


/3+ approaches —oo. This essentially means, the volume of the universe, i.e, 
does not hit singularity at any time. Similarly one can argue that the wave function 
dies out for x = 0 which does imply the probability of large /3_ is very small. But we 
note, Pr_|_ is time dependent and peak of Pr_|_ does shift to the right with time. Hence, 
the probability of having large /5+ is not small if we give long enough time for evolution. 
Hence, the volume can get arbitrarily large with time. 

3.2. Anisotropy 

The scalar found out of shear tensor can be written for a diagonal 

metric as, 

f 9u_ _ ^ ‘^ + ( — — —^ ^ 

\9'll 922/ \922 933/ \933 9ll / 

which yields 

= (30) 

for the metric (3) where are conjugate variables to /?±. 

4. Models with other Equations of State(0 < a < 1) 

In this section, we turn our attention to a 7 ^ 1. For a 7 ^ 1 we rewrite (15) with the 
following operator ordering, 

Jg^(l-a)(/3++4/3_)p^^g^(l-a)/3+^^^ _ gx/3(l-a)(/3++/3-)p^ gx/3(l-a)/3_^^ 

_48e-2^/3/^-e^T-“)(/5++2/3-) + 24pr] = 0. (31) 

The guiding principle behind choosing such operator ordering is to ascertain self- 
adjoint extension which is explained later in this section. 
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By replacing the momenta by the corresponding operators as described before, we 
have the Wheeler-DeWitt equation as 

_ ^(l-a)(/3++4/3_)_^ , x/3(l-a)(/?++/3-) ^ V3(l-a)0- ^ 

5/3+ 5/3+ dfi- dj3+ 

_48g-2x/3/3_gx/3(l-a)(/3++2/3_)]^ = 24z|^. (32) 

We now effect a transformation of variables as 


X+ = (i 

X- = e' 


^(l-a)/3+ 

v/3(l-a)/3- 


(33) 

(34) 


and use separability ansatz 4/ = 0(x+,x_)e to obtain 




1 1 d^cf) 

X- dxl ^ xl dx- 


48xh “x+V = 2AE(j). 


(35) 


This equation is not apparently separable in x± and thus one cannot investigate 
the behaviour analytically. But under a reasonable assumption, we can show the 
Hamiltonian given by (35) actually admits a self-adjoint extension. The assumption 
is that the Hamiltonian is a symmetric operator, by which we mean that the solution 
(if it admits any) obeys following conditions 



df 1" 
^X±Jo 


0 


(36) 


To facilitate a clear view, we can think of this assumptions in following fashion. 
Imagine, we obtain the solution to eq. (35) and restrict the solution space so as to satisfy 
eq. (36). Now what we call a reasonable assumption is that this restricted solution space 
is not null. This is actually quite natural in context of quantum mechanics, we assume 
similar conditions while we work on a free particle on a half-line or the whole real line. 
In case of whole line, 0 in the lower limit of the equation (36) is replaced by — cxd. 
The condition is reasonable in the present context, since, for example, any normalizable 
solution which vanishes at /3± = ±cx) will satisfy the condition. 

Now we can use Von Neumann’s theorem [28] asserting that a symmetric operator 
A dehned on domain V has equal dehciency index, if there exists a norm preserving 
anti-unitary conjugation map C \ V ^ T) such that [A, C] = 0, which, in turn, shows 
that A admits self-adjoint extension. 


It is easy to see that all the conditions for employing Neumann’s theorem are 
satished: 

(i) Complex conjugation map, C : "H) —)■ "H), between Hilbert space "H), takes </> to f*, 
which also belongs to the Hilbert space Ti). 


(ii) C is norm preserving , since whatever be the dehnition of norm, it involves fxp*, 
hence does not change under C. 
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(iii) C is anti-unitary and we have CHg = HgC. 

Hence, the theorem goes through and we have self-adjoint extension and thereby a 
unitary evolution for Kantowski-Sachs cosmology. The theorem can be understood in 
following manner as well. Since, satishes, 

Hgf+ = Z0+ (37) 

if and only if 

Hg(f)_ = —(38) 

where the map C induces a one-to-one map between two spaces, whose 

dimensions are actually named Deficiency Index and thereby make them equal, and 
we know that if a symmetric operator has equal deficiency indices, it does admit a self- 
adjoint extension. The detailed and rigorous proof can be found in [28]. It deserves 
mentioning that the extension may modify the boundary condition (36) by making it 
more strict (i.e the modified condition will imply eq. (36), not necessarily implied by 
eq. (36).) 

The role of operator ordering is crucial in the sense that with such an ordering we 
have the kinetic term multiplied with Hence, the condition (36) for its being 
symmetric is same as the condition for a standard Laplacian since the derivative with 
respect to x+ term is multiplied with y_ and vice versa. 


5. Discussion 

It has been shown that like anisotropic cosmological models with constant spatial curva¬ 
ture such as Bianchi I, V and IX, the Kantowski-Sachs cosmology can also be quantized 
where the problem of non-unitarity can be successfully eradicated. The trick is to hgure 
out the correct weight factor while dehning the norm and the inner product. This has 
been explicitly demonstrated for a stiff fluid given by P = p, i.e., for a = 1. We have 
calculated the shear scalar which is a positive definite quantity indicating that this 
unitarity restoration is not at the cost of anisotropy itself. In the case of a Bianchi I 
quantum cosmology, it had already been shown that the unitarity is not purchased at 
the expense of anisotropy [29]. For a 7 ^ 1, however, we could not solve the equation due 
to inseparability but we have shown that the Hamiltonian admits a self-adjoint exten¬ 
sion (with a suitable operator ordering). Thus the model can actually have a unitary 
evolution for any ideal barotropic (F = ap) fluid. 

One explicit example, namely that with a stiff fluid, along with the implicit proof 
for non-stiff fluids is good enough to disprove the folklore that anisotropic quantum 
cosmology necessarily involves non-unitarity. Here we have shown the Hamiltonian is 
indeed self-adjoint for a = 1, thereby has to admit a unitary evolution. Furthermore, 
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we have constructed explicit wave-packet showing unitary evolution explicitly. This 
restoration of unitarity is non-trivial in nature since the wave packet leads to a time 
dependent probability density, which integrates to a time independent hnite quantity, 
identihed as norm, leading to a conservation of probability. This scenario is in com¬ 
plete contrast with having energy eigenstates which trivially yields time independent 
probability density and thereby time independent norm irrespective of self-adjoint ness 
of underlying Hamiltonian. Energy eigenstates, thus, may not actually serve as a sig¬ 
nature of unitaritity. 

For a ^ 1, the proof is implicit in nature, it does not shed light on how to obtain 
a solution or how to do a self-adjoint extension. Yet it is good enough to consolidate 
the idea that non-unitarity is not generic to anisotropic models, it can be cured with 
suitable operator ordering or introducing a correct weight factor. 

As a bonus, it is also interesting to note that the stiff fluid model, a bonahde model 
which does not violate the classical energy condition, can actually give rise to a singu¬ 
larity free cosmology at the quantum level. This stems from the fact that scale factors 
have some intrinsic quantum fluctuation around its average value and the wavefunction 
goes to zero whenever scale-factor hits 0. 
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